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Section  1 
INTRODUCTION 


In  this  paper,  a  technique  is  developed  which  can 
be  used  to  estimate  the  average  reflectivity  of  ocean  bottom 
sediments  from  the  layering  statistics  of  the  sediment 
column.  Gilbert1  has  addressed  this  problem  with  Monte  Carlo 
simulation  techniques.  He  published  scatter  diagrams  of  the 
distribution  of  the  complex  reflection  coefficient  from  30 
realizations  of  a  randomly  layered  sediment,  showing  the 
reflection  coefficient  becomes  more  random  with  higher 
frequency.  Halthusen  and  Vidmar2  published  reflection  loss 
curves  for  several  realizations  of  randomly  layered  sediment. 
They  show  that  reflection  loss  can  be  greatly  reduced  by  the 
presence  of  layering.  In  this  paper  an  analytical  technique 
is  developed  for  calculating  the  probability  distribution  of 
the  complex  reflection  coefficient  for  a  random  ensemble  of 
layered-structures.  The  average  reflection  loss  calculated 
from  this  distribution  is  shown  to  agree  quantitatively  and 
qualitatively  with  observed  data. 

The  reflectivity  of  the  ocean  bottom  has  been 
measured  extensively  on  a  worldwide  basis  by  the  Naval  Air 
Defense  Center  and  the  Naval  Oceanographic  Office.  A  meas¬ 
urement  is  carried  out  by  detonatinq  SUS  (Sound,  Underwater 
Source)  charges  at  various  ranges  from  a  stationary  sonobuoy. 
The  energy  received  in  a  given  band  (typically  1/3  or  1 
octave  between  50  and  1600  Hz)  is  integrated  in  time  to 
obtain  the  total  received  intensity,  I.  The  bottom  loss  (BL) 
is : 

BL  =  -10  log  j Q( I)  +  SL  -  TL, 
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where  SL  is  the  source  level  in  dB,  and  TL  is  the  trans¬ 
mission  loss  in  dB  along  the  path  of  propagation. 

Near  normal  incidence  each  bottom-bounce  arrival  can  be  iden¬ 
tified  with  a  reflection  from  the  sediment,  no  refracting 
paths  being  present.  In  practice,  the  reflection  strengths 
are  often  quite  strong,  especially  when  the  sediments  are 
terrigenous  in  origin,  as  in  the  deep  ocean  abyssal  plains. 
The  observed  bottom  losses  generally  lie  between  3  and  10  dB, 
corresponding  to  pressure  reflection  coefficients  of  mag¬ 
nitude  .71  and  .32  respectively.  This  is  too  strong  to  be 
accounted  for  by  reflection  from  the  water-sediment  interface 
alone,  which  has  a  pressure  reflection  coefficient,  due  to 
density  and  sound  speed  discontinuities,  of  about  .2,  cor¬ 
responding  to  a  bottom  loss  of  14  dB.  In  this  paper,  it 
will  be  demonstrated  that  multiple  reflection  from  many  thin 
(~.25  -  1.  m)  sediment  layers  within  the  volume  of  the  sedi¬ 
ment  can  account  for  the  strong  reflectivity  of  terrigenous 
sediments . 


Section  2 


THE  BOTTOM  LOSS  UPGRADE  AND  THE  EXTRAPOLATION 
OF  GEOACOUSTIC  PARAMETERS 

The  Bottom  Loss  Upgrade3'4'5  is  a  database  of 
sediment  thickness  and  eight  parameters  which  characterize 
the  acoustic  properties  of  ocean  bottom  sediments.  The  data 
base  is  used  to  extend  the  sound  velocity  and  attenuation 
profiles  from  the  water  column  into  the  sediment,  as  shown  in 
Figure  1.  The  sound  speed  profile  depends  on  three  param¬ 
eters,  which  were  chosen  to  be  consistent  with  Hamilton's6 
canonical  profiles  for  terrigenous,  calcareous,  and  siliceous 
deep  sea  sediments.  The  attenuation  profile  is  a  linear 
function  of  depth.  The  two  parameters  of  the  attenuation 
profile  are  estimated  from  the  strength  of  bottom  refracted 
acoustic  arrivals.  Attenuation  values  agree  with  previously 
published  data  by  Mitchell  and  Fock.7  It  was  found  that 
attenuation  profiles  were  well  correlated  with  physiographic 
province . 


The  last  three  parameters  of  the  database  charac¬ 
terize  the  reflectivity,  or  equivalently  the  input  impedance, 
of  the  sediment,  as  described  in  Brekhovskikh . 8  It  was  found 
that  the  reflectivity  of  the  bottom  was  generally  stronger 
than  would  be  expected  from  the  density-sound-speed  discon¬ 
tinuity  at  the  water-sediment  interface.  The  strongest 
reflectivity  was  observed  in  abyssal  plains.  Figures  2  and  3 
are  histograms  of  the  minimum  and  maximum  impedances  of  the 
sediment  required  to  produce  reflection  strengths  of  the 
bottom  observed  in  the  NADC  data.  Figure  2  is  compiled  from 
data  sets  in  terrigenous  continental  rise  environments. 
Figure  3  is  compiled  from  data  sets  in  terrigenous  abyssal 
plain  environments.  If  the  sediment  is  treated  as  a  homo¬ 
geneous  fluid,  its  impedance  would  be  the  product  of  density 
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Figure  1.  Geo-acoustic  Parameters  of  the  Bottom-Loss 
Upgrade  Model 
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Figure  2.  Effective  Acoustic  Impedance  of  the  Ocean 
Bottom  for  Continental  Rise  Environments 
in  the  BLUG  Database 
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Figures  4a-4d.  Measured  Bottom  Loss  at  Four  Abyssal 

Plain  Sites  (pp.  3-12  and  3-13) 
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Section  5 

MODEL-DATA  COMPARISONS 


5.1  PREDICTIONS  AND  COMPARISON  WITH  DATA 

In  abyssal  plain  environments,  it  is  quite  common 
for  the  reflection  loss  to  be  in  the  5-  to  10-dB  range  for 
large  grazing  angle.  What  is  perhaps  more  surprising  is  that 
the  reflection  loss  often  decreases  with  frequency.  This 
phenomenon  is  referred  to  as  the  negative  frequency  depend¬ 
ence  of  bottom  loss.  Figures  4-a  through  4-d  are  examples  of 
bottom  loss  which  were  measured  in  four  abyssal  plains  at 
widely  separated  locations  in  the  North  Atlantic.  Above 
about  35°  grazing  angle  in  these  curves,  the  measured  bottom 
loss  becomes  fairly  level.  In  this  range  of  angles  the 
bottom  interacting  signal  consists  only  of  bottom  reflecting 
paths,  no  refracting  paths  being  present  at  short  range.  In 
each  case,  the  reflection  loss  is  much  less  than  the  1 5-dB 
level  attributed  to  reflection  from  the  water-sediment  inter¬ 
face.  Furthermore,  the  loss  at  1600  Hz  is  substantially  less 
than  the  loss  at  100  Hz.  This  is  the  negative  frequency 
dependence  phenomenon.  The  loss  curve  at  400  Hz  is  generally 
between  the  100-  and  1600-Hz  curves,  being  closer  to  the 
100-Hz  curve  at  the  lower  angles  and  falling  to  the  level  of 
the  1600-Hz  curve  at  the  higher  angles. 

A  prediction  of  average  bottom  loss  values  was 
carried  out  using  the  random  periodic  structure  algorithm. 
The  prediction  was  based  on  the  geoacoustic  model  in  Table 
1  : 


For  the  purposes  of  the  numerical  implementation,  a  discrete 
sine  transform  was  used  for  the  Uj  transform,  giving  an 
expression  for  the  pj  in  the  form: 

m  n  n 

p,(u.,9.)=  l  [  l  a.  sin  (lu  .  )cos  (k9.  )+  J  b  sin (lu .  )sin  (k 9 .  ) 1 . 
3  3  3  1=1  k=0  3  3  k=1  3  3 

(15) 

The  use  of  the  sine  transform  guarantees  the  boundary  condi¬ 
tion  pj(0,9j)  =  0  on  the  boundary  of  the  unit  disk. 


The  2nm+1  coefficients,  a.  .  and  b,  .  , 

kl  kl 

are  the  p^  in  Eq.  (9).  The  linear  relationship 
coefficients  of  pj^  and  p2, 


in  Eq.  ( 15) 
between  the 


P  2  ~  ^  i  P  i  t 


is  developed  by  evaluating  Eq.  (9)  at  2nm+1  discrete  values 
of  R2  =  e  i(u2+:‘-92^.  This  results  in  2nm+1  independent 
linear  equations  relating  the  p2  and  p^ 


The  mean  value  of  intensity  can  be  calculated  from 
p(u,9)  using  the  formula: 


<”  2n  . 

<I>  =  /  /  p(u,0)  e~4u  d  9  du  .  (16) 

0  0 
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periodic  random  structure  using  the  formalism  associated  with 
Eq.  (10). 


A  natural  coordinate  system  was  chosen  for  the 
implementation  of  Eq.  (8).  The  following  two  substitutions 
are  made,  for  the  reflection  coefficients  R2, 


R,  =  e-,U2+i62>  , 


and  for  RL  =  Lfc“ 1  (R2)  > 


Lt~l  ( R  2 )  =  e 


o-i  v  -(u  +i9  1 
-2lat  '  o  o  , 


where  u  and  0  are  chosen  so  that 
o  o 

-(uo  +  i0o)  R2-r 

e  =  - —  . 

1  -r  R2 

Then  Eq.  (8)  reduces  to 


p2(u2,02)  =  J0~1(R2)  /  p1(uo+Re(a)t,eo+Im(a)t) 

o 

-t/D 

exp(4 »Im( o)t  )  —  dt 


(13) 


In  this  coordinate  system,  the  probability  densities  can  be 
expanded  in  a  Fourier  sine  transform  in  uj  and  a  Fourier 
series  in  9j : 

®  n 

P-:  ( u  .  ,  9  .  )  =  /  sinU  u.)[  l  a.  (  5)cosfk0.  ) 

3  3  3  o  3  k=0  K  3 


n 

+  l  b  sin (k  0 .  ) ]  dC  .  (14) 


as  a  Fourier  Series,  Eq .  (8)  can  be  used  to  construct  a 

linear  map  from  the  expansion  coefficients  of  pL  to  those  of 

P2 : 


p2  =  Ai  Pi  ,  (10) 

where  At  is  a  matrix  and  p^  and  p2  are  the  expansion  coeffi¬ 
cients  of  px  and  p2  respectively. 

Similarly,  if  a  second  additional  layer  is  added,  a 
second  relation: 

P3  =  A2  P2  ,  (11) 

can  be  developed.  In  this  paper,  the  special  case  of  the 
periodic  random  structure  is  considered.  When  two  types  of 
layers,  each  with  its  own  acoustical  properties  and  thickness 
statistics,  are  alternated  then  the  coefficients  Pj  and  p3 
must  be  identical.  In  this  case,  Eqs.  (10)  and  (11)  reduce 
to: 


Pi  Aj  A[  Pj  .  (12) 

Eq.  (12)  was  implemented  numerically  by  iteratively  calcula¬ 
ting  the  eigenfunction  with  eigenvalue  closest  to  1.  When 
the  process  converged,  the  function  pt  corresponding  to  the 
expansion  coefficients  Pj  was  found  to  be  positive  definite, 
as  is  required  for  a  probability  distribution. 

An  actual  sequence  of  layers  in  ocean  sediment  is 
terminated  at  the  top  by  the  water-sediment  interface.  This 
final  terminating  layer  can  be  added  onto  the  top  of  the 
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choice  of  Q.  Thus  for  a  given  fixed  layer  thickness  t,  the 
following  relation  holds  between  pj  and  p2: 

P  2  (  R  2  » t )  =  p  !  ( L“  1  {  R  2  )  )  exp  (4  *Im(  a)  t  )  Jq-^Rj)  .  (7) 

The  random  layered  model  is  based  on  the  probabil¬ 
ity  distribution  of  the  reflection  coefficient  when  the 
layers  have  random  thicknesses.  This  can  be  obtained  by 
randomizing  the  density  in  Eq.  (7)  with  respect  to  the  thick¬ 
ness.  The  result  is  the  expected  value  of  p2  with  respect  to 
the  thickness: 


P  2  (  R  2  ^  —  /  P2(R2'^") 

0 


-t/D 


-t/D 


=  J0-1(R2)  /  Pi(l;1(R2))  exp  (4  *Im(  a)  t  )  - . dt, 

0 


where  the  expression  Lfcl(R2)  is 


L'‘(R2)  »  e 


-2 iat  R2-r 

*  ' 
1-r  R2 


and  D  is  the  mean  layer  thickness.  The  probability  density 
of  a  sequence  of  layers  of  random  thickness  can  be  calculated 
recursively  using  Eq.  (8),  by  adding  layers  to  the  top  of  the 
sequence  one  by  one. 

4.2  NUMERICAL  IMPLEMENTATION 

When  the  probability  densities  pi  and  p2  are 
expanded  in  terms;  of  a  set  of  basis  functions,  for  instance 


for  complex  reflection  coefficient  r,  where  r  is  the  complex 
conjugate  of  r.  Eq.  (2)  has  the  advantage  of  maintaining 
conservation  of  energy  and  reduces  to  Eq.  (1)  for  real  values 
of  r . 

When  an  additional  layer  of  specified  thickness  t 
is  added  to  the  top  of  the  sediment  sequence,  the  probability 
densities  for  Rx  and  R2,  p^Rj)  and  p2(R2rt)  respectively, 
are  defined  in  the  complex  unit  disk,  and  satisfy  the 
following  identity: 

/  p  2 ( R  2 , t )  dx  2  dy  2  =  /  P  i  ( R  i )  dx  t  dy x  ,  (3) 

Q  L-ifQ) 

where  Q  is  any  sufficiently  well-behaved  region  in  the  unit 
disk.  From  Eq.  (3)  it  is  possible  to  derive  an  algebraic 
relationship  between  pjtRj)  and  p2(R2,t).  By  simply  changing 
variables  in  the  second  integral,  there  results: 

/  p2(R2*t)  dx  2  dy  2  =  /  Pi  (L~l(R2))  JI1(R2)  dx  2  dy2  .  (4) 

Q  Q 

Here,  J~ 1  is  simply  the  Jacobian  of  1  ( R 2 )  which  is: 

,  0-kl 2) 

J"1  =  exp(4  »Im(  a)  t)  - 5 -  (5) 

1-r  | R 2 1  4 

=  exp(4*Im(a)  t)  J0-l(R2)  •  (6) 

It  is  a  simple  consequence  of  analysis  that  the  inteqrands  of 
Eq.  (4)  must  be  equal,  since  Eq.  (4)  holds  for  an  arbitrary 


coefficient  of  the  new  sequence  can  be  calculated  from  that 
of  the  previous  sequence  according  to  the  formula: 


f 


t* 


r2  - 


r  +  R^ 


1  +  rR  ,  e 


2i  at 

2T£t  =  Lt(Rl}  ' 


(D 


where 


Rt  is  the  reflection  coefficient  of  the  original 

layer  sequence, 

R2  is  the  reflection  coefficient  of  the  sequence 

with  the  addition  of  one  new  layer, 

r  is  the  reflection  coefficient  of  the  interface 
at  the  top  of  the  new  layer, 

t  is  the  thickness  of  the  layer,  and 

a  is  k  sin0,  the  wavenumber  times  the  sign  of  the 
grazing  angle. 

Attenuation  is  formally  added  to  acoustic  propaga¬ 
tion  through  the  use  of  complex  sound  speeds.  In  this  case, 
the  grazing  angle  0,  the  vertical  wave  number  a,  and  the 
reflection  coefficient  r  become  complex.  For  a  real  reflec¬ 
tion  coefficient  r,  Eq.  (1)  maintains  conservation  of  energy. 


that  is,  when  |Rt|  £  1  then  )  R  2 1  £  1.  However,  when  r 
becomes  complex,  Eq.  (1)  no  longer  conserves  energy,  and  can 
produce  reflection  coefficients  greater  than  one.  For  the 
purposes  of  this  paper,  Eq.  (1)  will  be  generalized  to: 


r  +  R,  e 


2  i  at 


R2  =  L(  R  i )  = 


1  +  r*Rje 
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Section  4 

DISTRIBUTION  OF  THE  REFLECTION  COEFFICIENT 


4.1  MATHEMATICAL  FORMULATION 

In  this  section,  the  probability  distribution  for 
the  reflection  coefficient  from  a  sequence  of  layers  of 
random  thickness  is  developed.  This  is  done  by  adding  layers 
one  by  one  onto  the  top  of  the  sequence.  At  each  step  of 
this  process,  the  distribution  of  the  reflection  coefficient 
Rn  for  the  n-layer  system  is  calculated  from  the  distribution 

of  the  reflection  coefficient  R  ^  of  the  n-1  layer  system. 
In  this  way,  the  acoustic  properties  of  fairly  complex 
layered  sequences  can  be  developed. 

In  this  paper,  the  case  of  the  periodic  random 
structure  is  analyzed.  In  this  case  the  sediment  consists  of 
two  materials  (i.e.,  mud  and  sand)  which  proceed  to  infinite 
depth  in  alternating  layers  of  random  thickness.  The  proba¬ 
bility  distribution  of  interest  is  the  distribution  of  the 
reflection  coefficient  over  the  entire  ensemble  of  possible 
layer  thicknesses.  The  advantage  of  the  periodic  random 
structure  is  that  if  one  additional  pair  of  mud  and  sand 
layers,  each  with  random  thickness,  is  added  to  the 
structure,  then  the  ensemble  of  possible  random  layered 
structures  constructed  is  identical  to  the  original  ensemble. 
Thus  the  statistical  properties  of  the  layered  system  are 
periodic,  even  though  some  individual  realization  of  the 
layered  structure  is  not. 


When  an  additional  reflecting  layer  is  added  to  the 
top  of  a  sequence  of  reflecting  layers,  the  reflection 


the  slump.  The  thickness  of  the  large  grain  layers  is  prob¬ 
ably  proportional  to  the  amount  of  sediment  on  the  conti¬ 
nental  shelf  available  for  transport,  which  also  increases 
proportionally  with  time. 

In  the  model  of  sediment  thickness  used  here,  the 
waiting  time  between  slump  events  is  assumed  to  be  a  random 
variable  with  exponential  distribution.  Accordingly,  the 
thicknesses  of  the  large  grain  and  small  grain  layers  are 
assumed  to  be  exponentially  distributed  random  variables. 
For  simplicity  it  is  assumed  that  the  thickness  of  each  layer 
is  independent  of  the  thickness  of  the  other  layers.  The 
exponential  distribution  function  has  the  explicit  form: 


F ( d )  =  1  -  exp{ -d/D) , 
with  the  associated  density: 

f (d)  =  ^  exp( -d/D) , 

where  d  is  the  random  variable  of  layer  thickness  and  D  is 
the  mean  layer  thickness.  The  large  grain  and  small  grain 
layers  are  assigned  different  mean  thickness  values,  the 
small  grain  layers  generally  being  the  thicker.  The  exponen¬ 
tial  distribution  is  often  referred  to  as  the  waiting  time 
distribution.  This  is  because  if  some  number  of  events  occur 
at  random  times  within  a  given  time  interval,  the  waiting 
time  between  the  individual  events  will  be  distributed 
according  to  the  exponential  distribution. 

For  the  purposes  of  calculation,  the  geo-acoustic 
properties,  such  as  sound  speed,  density,  and  attenuation, 
will  be  assigned  fixed  values  appropriate  to  the  type  of 
sediment.  While  the  actual  sound  speed  and  density  may  vary 
from,  say,  sand  layer  to  sand  layer,  the  variation  is  small. 
This  can  be  seen  in  measurements  published  by  Hamilton,17 
Horn,18  and  Tucholke.19 
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enough,  the  sediment  can  be  transported  to  great  distances. 
By  1965,  Potter  and  Schiedegger 1 6  demonstrated  a  good  corre¬ 
lation  between  the  thickness  of  the  graded  bed  and  the  grain 
size  in  turbidite  sediments.  They  showed  that  the  empirical 
relation  between  grain  size  and  bed  thickness  is  consistent 
with  predictions  of  the  rate  of  sediment  deoosition  based  on 
the  decay  of  turbulence  in  a  turbidity  current. 

3.2  THE  GEOPHYSICAL  MODEL 

At  frequencies  below  1600  Hz  sound  can  penetrate 
the  sea  floor  sediment  for  tens  or  even  hundreds  of  meters 
and  be  reflected  back  into  the  water  column.  Information  on 
sediment  layering  from  sediment  cores  is  generally  only 
available  for  at  most  about  10  meters.  Properties  of  sedi¬ 
ment  layering  below  the  core  can  at  best  be  characterized 
statistically.  In  this  section  a  simple  statistical  model  of 
sediment  layering  will  be  developed  based  on  the  sedimenta¬ 
tion  process. 

The  sediment  column  in  abyssal  plains  is  composed 
of  alternating  layers  of  large  and  small  grain  material.  The 
large  grain  material  is  laid  down  by  discrete  short-term 
events  associated  with  the  slump  of  sediment  off  the  conti¬ 
nental  shelf  and  transported  by  turbidity  currents.  Between 
the  slumping  events,  smaller  grain  material  is  laid  down 
continuously  by  the  more  gradual  processes  of  sedimentation 
in  the  deep  sea. 


In  both  types  of  sedimentation  process,  it  is 
reasonable  to  assume  that  the  thickness  of  the  layers  is 
proportional  to  the  time  between  the  slumping  events.  The 
finer  grain  layers  build  up  continuously  until  interrupted  by 


Section  3 

SEDIMENT  STRUCTURE 


3.1  NEAR  SURFACE  SEDIMENT  LAYERING 

Near  surface  sediment  layering  in  abyssal  plain 
sediments  has  been  observed  for  many  years  usinq  high  reso¬ 
lution  echo  sounders.  Sub-bottom  reflections  were  first 
observed  in  a  series  of  cruises  by  the  Woods  Hole  Oceano¬ 
graphic  Institution  between  1958  and  1962,  and  reported 
by  Hersey  and  Rutstein,9  Worzel,10  Heezen,  et  al.,11  and 
Menard.12  In  1959,  Worzel10  identified  a  particular  sub¬ 
bottom  reflection  in  the  southeastern  Pacific  with  an  ash 
layer  observed  in  sediment  cores. 

In  1965,  Hersey13  published  echo  sounding  from 
several  locations  in  the  North  Atlantic  Ocean  and  Mediter¬ 
ranean  Sea.  He  identified  observed  sub-bottom  reflectors 
with  alternating  layers  of  sorted  coarse  clastic  materials 
and  clays.  The  clastic  materials  were  either  sorted  sand  and 
silt,  laid  down  by  sediment  turbidity  currents,  or  volcanic 


According  to  Sadler,1"’  turbidity  currents  were 
first  proposed  as  a  mechanism  for  the  deposition  of  graded 
sediment  beds  by  Kuenen  and  Migliorini15  in  1950.  A  turbid¬ 
ity  current  is  a  suspension  of  sediment  in  seawater  which 
flows  down  the  continental  slope  and  spreads  out  on  the  ocean 
floor.  Settling  of  the  suspended  sediment  load  is  retarded 
by  turbulence  as  the  current  moves  down  the  slope.  Gradu¬ 
ally,  the  load  drops  out  in  an  abyssal  plain  as  the  current 
slackens.  If  the  initial  velocity  of  the  sediment  is  high 


(~1.5  g/cc)  and  sound  speed  (~1.5  km/sec),  approximately 
2.25.  As  can  be  seen  in  Figures  2  and  3,  the  required 
impedance  is  always  greater  than  2.25,  and  tends  to  be 
somewhat  higher  in  abyssal  plain  environments  than  in  the 
continental  rise  environments. 

During  the  development  of  the  Bottom  Loss  Upgrade, 
it  was  necessary  to  extrapolate  the  acoustic  properties  of  a 
given  region  of  the  ocean  bottom  to  other  regions  where  no 
data  was  available.  It  was  readily  apparent  from  the 
acoustic  data  that  abyssal  plains  exhibit  strong  reflec¬ 
tivity.  However,  the  actual  level  of  reflectivity  did  not 
seem  to  be  related  to  the  average  geophysical  properties  of 
the  sediment  column.  It  was  decided  to  extrapolate  the  geo¬ 
acoustic  parameters  from  area  to  area  based  on  the  observed 
or  surmised  presence  of  layering  in  the  sediment.  As  a 
primary  goal,  it  was  desirable  to  verify  that  this  was  a 
viable  technique.  The  optimum  solution  to  the  problem  would 
be  a  direct  method  of  estimating  acoustic  reflectivity  from 
the  observed  layering  in  available  core  data.  This  paper 
presents  preliminary  results  in  this  direction. 


Table  1 

TRIAL  GEOACOUSTIC  MODEL  FOR  BOTTOM  LOSS  PREDICTIONS 


MATERIAL 

DENSITY 

RATIO 

re/Water 

SOUND  SPEED 
RATIO 
re/Water 

MEAN  LAYER 
THICKNESS 

ATTENUATION 

dB/m/kHz 

Water 

1 

1 

N/A 

0 

Silty  Clay 

1  .5 

.996 

1 . 0  Meter 

.015 

Sand-Silt-Clay 

1  .7 

1  .004 

.25  Meter 

.015 

The  density  values  were  increased  sliqhtly  from  near-surface 
values  to  approximate  those  at  a  few  tens-of-meters  depth, 
although  bottom  loss  estimates  depend  only  on  the  ratio  of 
the  two  values. 

The  following  values  (in  Table  2)  of  mean  reflected 
intensity,  expressed  in  dB,  were  calculated  based  on  the 
values  in  Table  1 : 

Table  2 

ESTIMATED  BOTTOM  LOSS  IN  dB  BASED  ON  THE  GEOACOUSTIC 

MODEL  IN  TABLE  1 


O 

o 

Grazing 

Freq 

<Loss> 

m 

dB 

100 

7.49 

400 

7.68 

1600 

5.29 

80°  Grazing 

Freq  <Loss> 
Hz  dB 


.58 

.96 

4.64 


The  calculated  loss  levels  are  in  good  agreement  with  the 
four  cases  of  measured  loss  in  Figures  4-a  through  4-d.  At 
both  the  40®  grazing  and  80°  grazing  angles,  the  loss  at  1600 
Hz  lies  below  the  loss  at  100  Hz.  The  400-Hz  value  lies 
above  the  1600-Hz  value  at  the  40°  grazing  angle  and  moves 
down  to  the  1600-Hz  value  at  the  80°  grazing  angle.  For  this 
simple  case,  the  model  for  scattering  from  sub-bottom  reflec¬ 
tors  in  the  form  of  sediment  layers  of  random  thickness  seems 
to  account  for  both  the  observed  levels  and  the  neqative 
frequency  dependence  of  observed  bottom  loss  using  realistic 
values  of  geoacoustic  sediment  properties  and  layering 
statistics . 

5.2  LIMITATIONS 

The  algorithm  for  calculating  the  mean  value  of 
loss  based  on  Eq.  (8)  has  the  disadvantage  that  a  large 
number  of  terms  in  the  two-dimensional  series  expansion  of 
the  probability  density  must  be  calculated.  The  calculation 
of  the  N  coefficients  then  depends  on  iteration  using  an  NxN 
matrix.  The  large  memory  requirements  of  this  algorithm  may 
limit  the  number  of  terms  carried.  The  algorithm  has  yet  to 
be  extensively  tested.  In  particular,  the  sensitivity  of  the 
results  to  the  number  of  coefficients  is  unknown. 
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